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Abstract —In non-orthogonal multiple access (NOMA) down¬ 
link, multiple data flows are superimposed in the power do¬ 
main and user decoding is based on successive interference 
cancellation. NOMA’s performance highly depends on the power 
split among the data flows and the associated power allocation 
(PA) problem. In this letter, we study NOMA from a fairness 
standpoint and we investigate PA techniques that ensure fairness 
for the downlink users under i) instantaneous channel state in¬ 
formation (CSI) at the transmitter, and ii) average CSI. Although 
the formulated problems are non-convex, we have developed low- 
complexity polynomial algorithms that yield the optimal solution 
in both cases considered. 

Index Terms —5G, NOMA, fairness, outage probability, convex 
optimization. 

I. Introduction 

UTURE 5G communication systems need to support 
unprecedented requirements for the wireless access con¬ 
nection, targeting cell throughput capacities of 1000 x cur¬ 
rent 4G technology and roundtrip latency of about 1 msec. 
Towards this direction, three major 5G technologies named 
ultra-densification, millimeter wave, and massive multiple- 
input multiple-output, have attracted considerable attention in 
both industry and academia HI. In addition to these techno¬ 
logical implications, physical layer issues such as transmis¬ 
sion waveforms and multiple-access (MA) schemes should 
be reconsidered. A promising downlink MA scheme is the 
non-orthogonal multiple access (NOMA) which achieves high 
spectral efficiencies by combining superposition coding at the 
transmitter with successive interference cancellation (SIC) at 
the receivers El, 0. 

In contrast to the conventional orthogonal MA schemes 
(e.g. time-division multiple access (TDMA), etc.), NOMA 
simultaneously serves multiple users in the same degrees of 
freedom by splitting them in the power domain. Since NOMA 
is based on the SIC order, the served users achieve unequal 
rates and this could be critical for scenarios with strict fairness 
constraints. The work in a analyzes the performance of the 
NOMA scheme in terms of outage probability and achievable 
sum-rate for a fixed power allocation without discussing poten¬ 
tial fairness issues. In order to enable fairness, a introduces 
a cooperative phase after the NOMA downlink, where strong 
users relay data for the weak users but with the cost of 
extra channel resources (i.e., dedicated time slots are used 
for cooperation). On the other hand, fairness can be supported 
through appropriate power allocation (PA) of the superimposed 
transmitted data flows. 
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The main novelty of this work is the investigation of the im¬ 
pact of PA on the fairness performance of the NOMA scheme. 
In particular, we study the PA problem in NOMA from a 
fairness standpoint under two main system assumptions; i) 
when user’s data rates are adopted to the channel conditions 
(perfect channel state information (CSI)), and ii) when users 
have fixed targeted data rates under an average CSI. Although 
the resulting optimization problems are both non-convex, we 
develop low-complexity bisection-based iterative algorithms 
that provably yield globally optimal solutions. We further 
illustrate that each iteration subproblem can be optimally 
solved in closed and semi-closed form for the two cases, with 
no specialized optimization software. The results show that 
the NOMA scheme outperforms conventional MA approaches 
by significantly improving the performance of the worst user. 

II. System Model 

We assume a single-cell downlink topology consisting of 
one base station, B, and N users, Ui, with i G Af = 
{1,..., A^}; all terminals are equipped with a single antenna. 
B has always data to transmit for each user (saturated scenario) 
and its total available transmitted power is equal to P. All 
wireless links exhibit independent and identically distributed 
(i.i.d.) block Rayleigh fading and additive white Gaussian 
noise (AWGN). This means that the fading coefficients hi 
(for the B ^ Ui link) remain constant during one slot, but 
change independently from one slot to another according to a 
complex Gaussian distribution with zero mean and variance 
the variance captures path-loss and shadowing effects. 
Without loss of generality, the channels are sorted as 0 < 
1^1 P < < ■ • ■ i-e., the i-th user always holds the 

i-th weakest instantaneous channel. The AWGN is assumed to 
be normalized with zero mean and variance cr^. 

A. NOMA scheme 

The NOMA scheme allows B to simultaneously serve all 
users by using the entire bandwidth to transmit data via a 
superposition coding technique at the transmitter side and 
SIC techniques at the users lb) Eq. (6.25)]; in this case, user 
multiplexing is performed in the power domain. According to 
the NOMA principles, B transmits a linear superposition of 
N data flows by allocating a fraction Pi of the total power to 
the i-th data flow. Each receiver employs a SIC technique and 
is able to perfectly decode the signals of the weakest users 
i.e., the i-th user can decode the signal for the m-th user with 
m < i, removing completely inter-user interference. In this 
way, the )-th user can remove interference from the weakest 
users and its achievable rate associated with the decoding of 
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the m-th data flow is given by 141 
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where Rm,m < Ri,m, m < i due to the ordering of the channel 
coefficients. 


B. Outage probability performance for NOMA 

In case B does not have any instantaneous channel feed¬ 
back, it transmits with a target spectral efficiency tq bits per 
channel use (BPCU) for each data flow and an appropriate 
performance metric is the outage probability. An outage event 
occurs at the i-th user when is not able to decode its own 
data flow or the data flows of the weakest users m < i. By 
using high order statistics and following similar steps to a, 
the outage probability for the i-th user is 


^NOMA 
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where 7 ^,^ = -j-. ^ 

X{N - i + 1 +k), Ci = niax{Ci,C 2 , ■ • ■ ,CJ> ^ = 1 /ct^, 

-hP^> fo EL+i Pi (see 

Eq. ( 8 ) in 0), and (l 2 bl i follows from the binomial theorem. 
It is worth noting that here, we give a closed form expression 
for i.i.d. Rayleigh fading, while the analysis in a refers to a 
different channel probability distribution. 


III. Fairness for NOMA systems 

The NOMA scheme enables a more flexible management of 
the users’ achievable rates and is an efficient way to enhance 
user fairness. In this section, we consider two fundamental fair¬ 
ness criteria that refer to NOMA systems with instantaneous 
and average CSI, respectively. 


A. Max-Min fairness with instantaneous CSI 

If a continuous channel feedback is available at the trans¬ 
mitter side, users’ rates can be allocated according to their 
instantaneous channel conditions. In this case, a suitable 
criterion is the max-min fairness that maximizes the minimum 
achievable user rate and is formulated as follows 


max min 
(3 i£N 

(3a) 

N 


s.t. '22 Pj ^ 

(3b) 

7=1 


0 < Pj, for j G AC. 

(3c) 


Algorithm 1 ; Optimal solution to problem Q 

1: Init. tLB = 0, tjjB = log (^1 -f 
2 : while (tuB - tLB > f) do 

3: Set t = {tuB + tLB)/‘2', Solve LP ® to obtain f3^^. 

4: if I ^2 Pi^ — ^ 

5: Set tLB =t; p* = (3^^■, r* =t. 

6: else 

7: Set tijB = t. 


Problem 0 is not convex and hence hard to solve directly 
using standard optimization solvers. In this section, we trans¬ 
form the problem into a sequence of linear programs (LPs) and 
develop a customized low-complexity polynomial algorithm 
for its optimal solution. Towards this direction, we state 
proposition [T] 

Proposition 1. Problem 0 is quasi-concave. 


Proof: A maximization optimization problem is quasi¬ 
concave when the objective function is quasi-concave and 
the constraints are convex. Clearly, the constraints of Q are 
convex since (l3bl) and are linear. For the objective function 
to be quasi-concave, all its sublevel sets must be concave a, 
i.e., St = {mini?^°'’^^(/3) > t}, for f G M, which denotes the 

i ’ 

set of P for which the objective function is larger than t. Due to 
the min operator, it is true that St = >t,i€ A/"}. 

Nonetheless, set St is concave for f G R, as constraints 
7 jNOma(^) >t,i € ff can be expressed as 


P^P\h^\^ > (2‘ 


/ N 

P\h.\^ Pi 

V Z=i-|-1 


eA/", (4) 


which are linear inequalities and hence concave, completing 
the proof. ■ 

Let r* denote the optimal objective function value to quasi¬ 
concave problem Q. For a specific constant value t, if the 
LP 


Find P subject to constraints Obl i. d^ . 01 (5) 


is feasible then r* > t, otherwise r* < t. Equivalently, one 
can solve the following LP 


m\n2^Pi subject to constraints d^ and 0 - 

^ iGjV 


(6) 


and check if the solution satisfies (l3bl i. This implies that 
by appropriately bounding t through a bisection procedure 
(Algorithm [T]i, the optimal solution to a, within a desirable 
accuracy e, can be obtained by solving a sequence of feasibility 
LPs of the form (| 6 ]l. 

Although a is LP and can be solved with standard op¬ 
timization solvers, its solution can be obtained in closed- 
form with linear computational complexity from a customized 
algorithm. Towards this direction Proposition |2] is essential. 


Proposition 2. The optimal solution of satisfies all 
constraints d3^ with strict inequality (inactive constraints) and 
0 with equality (active constraints). 
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Proof: Because the problem is convex, the following 
Karush-Kuhn-Tucker (KKT) conditions are necessary and suf¬ 
ficient for optimality of Q, for i £ J\f: 


+ Afc(2‘ - l)P|/ife|2 + 1, 


k<i 


> (2‘-1) (pM^ Y. + 

\ 1=1+1 ) 


> 0, A, > 0, /r, > 0, 


A, 1^(2* - 1) (^P\h,\^ ^ ^ 


(7) 

( 8 ) 
(9) 


= 0 , 




( 10 ) 

( 11 ) 


where Ai and fj,i, i € Af, are the Lagrange multipliers for 
constraints (HI and (O respectively. The right hand side (r.h.s.) 
of ® is strictly positive for all i £ Af, as > 0, f > 0 and 
/3i > 0; hence, the left hand side (l.h.s.) has to be strictly 
positive which implies that /3i > 0 and Hi = 0 (due to (fTTI) l. 
which completes the first part of the proposition. In a similar 
fashion the r.h.s. of (|7]) is strictly positive, which implies that 
Ai > 0, i £ Af, as fii = 0. Since Xi > 0, condition (fTOl i 
implies that all constraints (O must be enforced with equality 
which completes the proof. ■ 

The intuition behind proposition |2] is that a\\ /Si, i € Af 
must be positive for the problem constraints to be satisfied, 
while all Ri^i must be equal at the optimal solution. Based on 
proposition |2] the optimal solution to LP (|6]l can be obtained 
from Theorem [T] 


Theorem 1. The optimal solution to I© is given by: 
‘yt — 1 f ^ A 

^‘p\h+Y,D>+^l 

V l=i+l 


P^ = 


P\h, 


= ( 12 ) 


Proof: The proof is a direct consequence of the active 
inequality constraints (IHi from Proposition |2] The analytical 
solution emanates from the fact that normalized power Pi 
only depends on the power allocated to the stronger channels 
i -f 1,..., iV. Hence, the optimal objective function value can 
be computed by allocating power from the strongest to the 
weakest channel in succession. ■ 

Notice that due to the special structure of the problem, the 
optimal solution can be obtained in closed form (see (fT^ l. 
with all users having data rate equal to t. In addition, it can be 
easily observed that the computational complexity of solving 
(|6]l is 0{N), i.e., linear to the number of users. 


B. Min-Max fairness with average CSI 

Optimizing the performance of communication systems 
under average CSI information is an important and challenging 
problem with practical interest. Towards this direction, we 
propose to optimize the outage probability of all users under 
knowledge of the distribution and order of channels in NOMA 
systems. The considered optimization problem is defined as 

minmax (13a) 

/3 ieAT 

subject to constraints Obl i and (O, (13b) 


where 7 :>noma^^^ given by (l2bT i. Naturally, increasing 
the total power monotonically decreases the optimal outage 
probability of ( fTSl l. Hence, a bisection procedure similar to 
Algorithm [T] can be followed, where each step involves the 
solution of the following problem 

min^/3i, (14a) 

s.t. <t, Pi> 0. (14b) 

For a particular value t, if the solution of Gil satisfies 
Szga/" Pi — ^ '^hen the optimal solution of is r* < t, 
otherwise r* > t. Hence, the optimal solution of Gl . can be 
obtained by optimally solving a sequence of Gl . Nonetheless, 
problem G1 is non-convex due to the presence of terms 
exp{—5i^kCi) in G4al i which are neither convex or concave. 
To simplify the problem, the next proposition allows the 
elimination of 

Proposition 3. At the optimal solution of problem G4b . it is 
true that \hip = C,i < ... < \hfq\^ = Cat. 

Proof: The proposition will be proved by induction. 
Assume that the optimal solution of G1 is /3*. From the 
definition of NOMA outage probability (l2at for i = 1, 
we have that P{|/ii|2 > Ci} =P{|/rip > 

Since /3* only affects the outage probability of user 1, we 
want to select Pi as small as possible to minimize the total 
used power which implies that P{|/iip = Ci}. Assume that 
\hip = (i,...,\hkp = Ck', we will prove the result for 
fc + 1. From (I2al) we have that ]P{nm=i.....fe-i-i > Cm} = 

P {\hk+i\'^ > Cfe+i}> as from the above assumption is it true 
that \hk+ip > Cm = m = Because Pk+i 

does not affect the outage probability of users A: + 2 to A^ it 
can be selected to minimize the total consumed power which 
is true when \hk+ip = Cfc+i- This completes the proof. ■ 
The intuition behind proposition [3 is that the target spectral 
efficiency needs to be satisfied with equality for all users. 
Proposition [3 implies that at the optimal solution of G4t . it 
is true that Ci = Ci- Next, we show that despite of GU being 
non-convex, its globally optimal solution can be obtained in 
semi-closed form. 

Theorem 2. The optimal solution of Gil is given by: 

2 - N 

P: = ^+fo Y PkA = N,...,f (15) 

k=i+l 

where Q, i G Af denotes the solution of the 1-D equation 

7z,fc (1 - exp(-5,,fcCz) ] = (16) 

Proof: To prove the above theorem, variable transforma¬ 
tion Ci = ^ is made on problem G1 yielding 

min Y o'7o(l + /{PC,i), (17) 

ieN 

Ci>0, i €J\f. 
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Fig. 1. Optimal NOMA fairness rate for different Fig. 2. Comparative NOMA and TDMA fairness Fig. 3. Comparative outage probability per- 
number of users and total power under instanta- rate perfoimance under instantaneous CSI for formance between NOMA, TDMA and fixed 
neous CSI. different number of users and total power. NOMA under average CSI. 


Interestingly, the constraints in problem (fTTl i are uncoupled 
and the objective is additive; hence, the problem can be 
decomposed into N one-dimensional problems, where the z-th 
subproblem is of the form 


-crlfoil + roT V(-PCi), 

(18a) 

li,k ( 1 - exp(-(5i,feCi)) < 

(18b) 



c. >0. 

(18c) 


Although, subproblems ( fTSl l are non-convex (due to dlSbl l). it 
can easily be observed that dlSab is monotonically decreasing 
in Q, while constraint dl8bb is monotonically increasing in 
(i. Hence, the optimal solution to (fTSb is the largest value of 
Q that satisfies dl8bb . which is obtained when the constraint 
is active. Having obtained the optimal solution of each sub¬ 
problem the optimal solution to (3* can be obtained by 
back-substitution yielding ( fTSb . completing the proof. ■ 

The solution of each (fTbl l can easily be performed using 
Newton’s or bisection method in 0{Nlog{e)), where e is the 
required accuracy. Once C* has been obtained, /3* is computed 
from (flSb in 0{N), yielding a total computational complexity 
of 0{N‘^log{e)) or 0{Nlog{e)) if solved in parallel. Note 
that the methodology proposed to optimally solve (fT4b . can 
be applied to any underlying outage probability distribution. 

IV. Numerical Results 

We evaluate the performance of the developed algorithms 
by solving 1000 randomly generated problems for different 
parameter configurations. All problem instances follow the 
system model introduced in Section |II] with — 1 and 
= 1. The fairness performance of the conventional TDMA 
scheme is used as a benchmark since it refers to the orthogonal 
allocation of the available degrees of freedom and is equivalent 
to any orthogonal multiple access scheme ( El, Sec. 6.1.3); 
for the instantaneous CSI case, the optimal TDMA allocation 
can be solved by using a methodology similar to Algorithm 
[T] while the optimal outage probability is obtained for equal 
time-sharing and power-split. 

Fig. 1 demonstrates the achievable maximum fairness rate 
for different P and N. As expected, increasing P or reducing 
N improves the achievable fairness rate. Interestingly, the 


performance gain from N = lOtoW = 5 is significantly 
higher than from W = 20 to = 10. Notice also that as P 
increases the rate of improvement reduces because the fairness 
data rate is a logarithmic function of power. 

Fig. 2 provides a comparison between NOMA and TDMA 
by depicting the achievable fairness rate for different P for 
varying number of users. Clearly, NOMA is significantly better 
than TDMA. In fact, as the number of users increases, the 
advantage of NOMA over TDMA increases almost linearly. 
NOMA is also substantially better than TDMA in terms of 
computational complexity, as TDMA requires the solution of 
a sequence of convex programs. 

Figs. 3 compares the outage probability between NOMA 
and TDMA under average CSI information as a function of the 
total transmitted power for target spectral efficiency tq = 0.05 
BPCU and rg = 0.50 BPCU, for W = 5. In TDMA the 
optimal outage probability is obtained for equal time-split ratio 
and power among users. Results are also shown for a fixed 
NOMA PA scheme proposed in iHJ, with Pm = 
m G M, where p is selected such that 
can be easily observed that in all cases considered, NOMA 
outperforms TDMA by an order of magnitude. NOMA also 
has at least five times better performance compared to the 
fixed NOMA PA scheme. As expected, higher target spectral 
efficiency results in worse outage probability because it is 
more difficult to be satisfied. For validation purposes, outage 
probabilities were obtained numerically for 10® problems 
instances illustrating maximum absolute relative error 0.7% 
for all cases considered. 

V. Conclusions 

In this paper, the problem of optimal PA to maximize fair¬ 
ness among users of a NOMA downlink system is investigated 
in terms of data-rate under full CSI and outage probability 
under average CSI. Although the resulting problems are non- 
convex, simple low-complexity algorithms are developed that 
provide the optimal solution. Simulation results demonstrate 
the efficiency of NOMA, achieving fairness performance that 
is approximately an order of magnitude better than TDMA 
in the considered configurations. The main results of this 
work show that NOMA can ensure high fairness requirements 
through appropriate PA and is a promising MA scheme for 
future 5G communication systems. 
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